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1.  INTRODUCTION 


When  small  panicles  suspended  in  a  gaseous  medium  flow  through  a  long  cylindrical  tube  or 
conduit,  the  random  Brownian  movement,  or  diffusion  of  the  particles  which  may  be  submicron 
aerosols,  atoms,  or  ions,  may  bring  them  into  contact  with  the  walls,  where  they  adhere  or  lose  their 
charge.  From  the  fraction  of  particles  penetrating  the  conduit,  the  diffusion  coefficient  and  particle 
size  may  be  calculated  with  the  use  of  the  appropriate  relations.  The  case  in  hand  is  when  no  particles 
enter  the  conduit,  and  formation  in  flight  occurs  within  the  containment;  this  arises  when  air 
containing  a  radioactive,  rare  gas  enters  the  conduit  through  a  high-efficiency  filter  located  far 
upstream  from  the  tube  configuration.  As  the  radioactive  gas  flows  through  the  geometry,  it  decays, 
giving  rise  to  the  steady  production  of  a  certain  number  of  daughter  atoms  per  unit  volume.  Unlike 
the  radioactive  gas,  the  daughter  atoms  adhere  to  the  wall  and  are  lost  by  diffusion.  This  dispersion  of 
atoms  may  be  considered  to  be  an  atomic  aerosol  (Tan  and  Hsu  1971)  since  it  has  the  same  property 
as  submicron  aerosols  in  that  they  can  be  collected  at  a  surface.  The  radioelements  diffuse  to  the  tube 
walls  where  they  decay  into  other  radioelements. 

Formation  in  flight  diffusion  equations  have  been  derived  for  cylindrical  tubes  (Tan  and  Hsu 
1971;  Berezhnoi  and  Kirichenko  1964)  and  for  flat  channels  (Berezhnoi  and  Kirichenko  1964), 
including  the  mass  transfer  of  aerosols  with  axial  diffusion  in  narrow  rectangular  channels  (Tan  and 
Hsu  1972). 

Swirling  flows  related  to  ducts  have  been  studied  extensively  (Fromm  1963;  Pao  1967;  Textor 
1968).  A  well-defined  configuration  of  a  confined  axially  decaying  vortex  flow  was  introduced  in  a 
study  by  Lavan,  Nielsen,  and  Fejer  (1969).  Also,  Tung  and  Soo  (1973)  presented  a  study  applying  the 
fluid  phase  relations  of  Lavan,  Nielsen,  and  Fejer  (1969)  to  the  case  of  a  gas-solid  suspension. 

The  present  study  utilizes  the  geometry  of  Lavan,  Nielsen,  and  Fejer  (1969)  (see  Figure  1)  to 
study  the  case  of  a  laminar  swirling  flow  of  atomic  particles  formed  in  flight  under  electrostatic  field 
effects.  The  fluid  phase  is  assumed  incompressible  and  fully  developed  at  both  far  upstream  and 
downstream  positions  from  die  juncture  of  the  two  pipes.  A  forced  vortex  is  generated  by  the  rotating 
[ripe,  and  it  decays  as  the  flow  passes  through  the  stationary  duct  because  of  wall  surface  friction.  At 
far  downstream,  the  swirl  vanishes,  and  it  is  assumed  that  the  axial  velocity  converts  back  to  the  fully 
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developed  laminar  parabolic  profile.  The  panicle  concentration  is  assumed  low  enough  so  that  the 
effect  on  the  fluid  phase  due  to  the  presence  of  the  panicles  is  negligible  (Soo  1967).  The  solutions  of 
the  fluid  phase  velocities  can  then  be  effectively  utilized  to  solve  for  the  paniculate  relations.  The 
fluid  phase  velocities  are  obtained  by  solving  the  Navier-Stokes  equations  numerically. 

2.  FORMULATION 

The  motion  of  the  fluid  phase  is  given  by  solving  the  continuity  and  the  Navier-Stokes  equations 
for  a  steady,  laminar,  incompressible,  and  axisymmetric  swilling  flow.  The  solution  is  given  in  terms 
of  the  radial,  tangential,  and  axial  components  of  the  fluid  velocity  (u,  v,  w  in  r,  0,  and  2  coordinates, 
respectively),  pressure  p,  the  density  p,  and  v,  the  kinematic  viscosity  of  the  fluid  phase  material. 
Hence,  the  fluid  phase  relations  arc  (Schlichting  1968): 


and  the  continuity  relation  is 


(ru)  +  ^  (rw)  -  0. 


(4) 


The  expression  describing  the  steady-state  mass  diffusion  of  a  constituent  in  a  generating  but 
nonreacting  binary  gas  mixture  under  electrostatic  influences  flowing  through  the  system  described  in 
Figure  1,  assuming  azimuthal  symmetry  and  constant  coefficient  of  diffusion  is 


where  F  is  the  inverse  of  relaxation  time  constant  defined  as 

F  *  9jl  /  2a1  c  ,  (7) 

where  ji  and  c  are  the  viscosity  of  the  fluid  and  material  density  of  the  particulate  phase,  respectively, 
and  a  is  the  radius  of  a  particle. 


The  boundary  conditions  of  the  solution  take  the  following  forms: 


.  r,  dw  dc  d<b  „ 

atr«0:«»v  -  «  __  ■  »  0 ; 

dr  or  dr 

(centerline) 


at  r  -  Rg :  u  -  w  ■  c  »  $  ■  0;  v 
(wall) 


RgG,  r$0 
0,  0  <  z  ^  +  <*> 


(9) 
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at  z  -  -oo;  u  -  c  -  $ 
(far  upstream) 


0;  v  ■  rft,  w 


(10) 


at  z  ■  +oo;  u  «  v  ■  0;  w 
(far  downstream) 


(ID 


where  ft  is  the  constant  angular  velocity  of  the  rotating  pipe,  w  is  the  mean  axial  flow  velocity,  and  Ra 
is  the  radius  of  the  pipe.  • 


The  flow  configuration,  together  with  the  boundary  conditions,  is  shown  in  Figure  1.  The  pipe 
section  at  the  left  of  the  r-axis  rotates  with  constant  angular  velocity  ft  and  creates  the  swirl,  while  the 
pipe  section  at  the  right  is  held  stationary.  At  z  =  0,  the  two  pipes  join  smoothly  with  no  effect  on  the 
flow.  At  far  upstream  (z  =  -«)  and  downstream  (z  =  +«>),  the  flow  is  assumed  fully  developed. 


Since  the  radioelements  are  completely  annihilated  at  the  walls,  c  *  0  at  r  *  ft,;  also,  because  of 
the  grounded  walls,  $  =  0  at  r  =  /?„  Furthermore,  it  is  assumed  that  c  =  0  at  z  = 


The  following  nondimensional  forms  may  be  introduced  as: 

R  -  r/R0,  Z  -  z!R0, 


W  -  JL,  U  -  JL,  V  -  v/(ft  ft); 
w/2  w/2 


C  F-p/(pft,2ftJ); 

C*  -  Dpc/(oft?),  4>*  -  q  */[mpFDp) ; 
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where  is  the  stream  function  at  the  wall.  The  Reynolds  number  is  given  by 


Rey  * 


2 wRa 

O 


V 


(12) 


and  the  swill  ratio  is  given  by  the  ratio  of  the  tangential  velocity  of  the  rotating  pipe  to  the  mean  axid 
velocity, 


S  * 


(13) 


The  diffusive  Peclet  number  is  the  ratio  of  the  inertial  to  diffusive  forces, 


p 

•  D, 


and  the  electrostatic  charge  parameter  is  given  by 


1 

(  v 
JL 

4 

m . 

v  p) 

e„FD 

O  p 


(14) 


(15) 


The  classical  numerical  approach  which  converts  the  velocity  components  into  stream  function  f, 
vorticity  £,  and  circulation  T  was  chosen  because  of  its  well-established  stability. 

The  relations  for  the  fluid  phase  after  they  have  been  expressed  in  nondimensional  form  were  cut 
into  explicit  finite  difference  molecules  (Southwell  1940;  Ames  1972)  and  solved  by  numerical 
relaxation  yielding  results  consistent  with  Lavan,  Nielsen,  and  Fejer  (1969)  and  Tung  and  Soo  (1973). 
As  a  result  of  the  axial  symmetry  of  the  flow  field,  only  the  flow  in  the  rectangular  region  defined  by 
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D  ■  [(/?,X)  |  0£  1  andOSXS  l]  <16> 

need  be  considered  where  X  is  a  mapping  function  given  by  Tung  and  Soo  (1973)  and  used  to 
transform  the  z-axis  into  the  x-axis.  On  D  itself,  parallel  mesh  points  (41  x  41)  were  uniformly  spaced 
by  an  amount  A  in  both  radial  and  axial  directions.  Initial  assumed  values  for  T,  £,  and  ¥  were 
specified,  then  an  iterative  procedure  consisting  of  sweeps  of  the  interior  mesh  points  was 
implemented  until  a  convergence  criteria  (Forsythe  and  Wasow  1967)  was  satisfied. 

The  diffusion  and  Poisson  equation  can  be  written  in  nondimensional  form  as 


(  > 

r. 

( 

I2  J 

i 

V 

2/^*  i/’*  \ 

+  X’  w 


j_  ac*  ^  a2c*  ^  a<i>*  a c* 
r  +  a/?2  +TF"5F 


-4aC*2 


+ 1, 


(17) 


i8**  + 

R  "3F+  d  R1 


-4a  C* 


(18) 


subject  to 


C*  (R,  X«0)  -  9*  (/?,  X-0)  «0; 


(19) 


C*  (1,X)  -  9*(1,X)  -  0; 


(20) 


dC*  n  (21) 

(0,X)  *  (0, X )  -  0 . 

Defining  F(X)  as  the  ratio  of  the  total  particle  flux  over  a  cross-section  at  distance  X  from  the  far 
upstream  position  (at  z  =  -«>)  to  the  rate  of  formation  of  the  radioelements  in  the  same  element  of  the 
tube,  we  have 


1 

F(X)  -1  f  WC*RdR. 

x  6 


(22) 
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After  the  fluid  phase  is  solved,  the  respective  velocities  can  be  implemented  into  the  solution  of 
the  differential  system  encompassing  Equations  17  to  21.  Due  to  the  universal  stability  of  the  implicit 
method,  this  parabolic  system  was  expressed  in  terms  of  the  implicit  finite  difference  discretization 
technique  such  that  the  diffusion  and  Poisson  equations  take  the  forms,  respectively, 


<»(jc  i+u**hjc 


(23) 


yl .  .A*  .  ,  +  y2.  .A*  .  +  v3  .  .A*  .  ,  +  4a C  .  , 
1  i,j*  j-i  1  i,j *  j  1  i,j*  y+i  i+i 


o. 


(24) 


where 


\ 

j 


i  +  1, 

2AR 


+ 


1 

2RjAR 


1  . 
(A  R)2  ’ 


(25) 


V2  y 


AX  (AR)2  1  'J 


(26) 


Vi,j 


(  \ 


\2  J 


U 


i  +  l  ,j  1.0 


2  A  R  2Rj(AR)  (A  R) 


2  ' 


(27) 


1 

Far 


f  c*  c * 

c  /  j+i  -  0  f,/-i 


./  ~  2AR 


TEW 


fr*  C *  N 

C  / ,)+ 1  -  C  i  J- 1 


(28) 


2AR 


(29) 


/ 


8 


(30) 


ai 


i  J 


(  \ 

p 


X'i+lWUl,JC  ij  +  ,  0 
AX 


yi 


* "  7*W  "  ^JaxWj  ; 


y2i’Ja~7*R?; 


Y3/  i  - 


lf  *j  "  Jar?  +  1*ji**)  ‘ 


(31) 

(32) 

(33) 


Convergence  (Forsythe  and  Wasow  1967;  Hombeck  1975)  was  satisfied  for  all  Peclet  and 
electrostatic  charge  parameters,  with  respect  to  the  mesh  sizes  considered  (AX,  A R  =  0.01, 0.025, 0.05). 
The  sparse  and  unsymmetric  system  resulting  from  the  finite  difference  expressions  for  the  diffusion 
and  Poisson  system  was  solved  by  Gaussian  elimination  with  full  and  partial  pivoting  with  the  aid  of 
the  Crout  reduction  technique. 


3.  RESULTS  AND  DISCUSSIONS 

Since  the  fluid  phase  relations  have  already  been  solved  (Lavan,  Nielsen,  and  Fejer  1969;  Tung 
and  Soo  1973),  they  will  not  be  discussed,  except  for  the  precision  with  which  they  matched  the 
available  comparable  calculations,  which  was  within  0.01%  for  most  cases. 

The  effect  of  the  electrostatic  field  parameter,  a,  on  the  fraction  of  penetration  over  the  complete 
axial  distance  from  far  upstream  to  far  downstream  positions  can  be  seen  in  Figure  2.  Clearly,  the 
increase  in  charge  causes  a  decrease  in  penetration,  F(X),  since  more  atomic  particles  are  attracted  to 
the  walls  where  they  become  completely  annihilated. 
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Figure  3  represents  the  effect  cf  the  charge  on  the  centerline  concentration  over  the  complete 
axial  distance.  A  decrease  in  profile  is  accompanied  by  an  increase  in  charge  since  the  particles  are 
lost  at  the  boundary.  Figure  4  identifies  tire  greater  electrostatic  potential  due  to  the  increase  in 
electric  field. 

Hence,  the  electrostatic  charge  parameter, «,  has  a  significant  impact  upon  nuclear  particle 
deposition  in  vortex  flow  when  1  £  Rey  £  10, 1  £  S  £  12, 1  £  P,  £  10,  and  1  £  <*  £  10,  which 
represent  an  appreciable  realistic  range.  Increasing  «  increased  deposition  and  potential,  regardless  of 
Rey,  P.  and  S  variations  within  these  ranges. 


10 


o 


BO 


N0LLVHIN30N03  3NriH3IN3D 


2 


Rgurc  3.  AxW  Distribution  of  the  Cnneriine  Concentration  at  Different  Electrostatic  Fields. 
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LIST  OF  SYMBOLS 


a 

c 

C * 

*>, 

F 

m, 

P 

P 

P. 

Q 

r 

R 

Rey 

R0 

S 

u 

U 

v 

V 

w 

W 

W 

X 


-  Radius  of  a  particle 

-  Particle  concentration 

-  Dimensionless  particle  concentration 

-  Particle  diffusivity 

-  Inverse  of  relaxation  time 

-  Mass  of  a  particle 

•  Fluid  static  pressure 

-  Dimensionless  pressure 

-  Peclet  number 

-  Rate  of  particle  generation  per  unit  of  volume 

-  Radial  coordinate 

-  Dimensionless  radial  coordinate,  r/R, 

-  Reynolds  number 

-  Pipe  radius 

-  Swirl  ratio 

-  Radial  velocity 

-  Dimensionless  racial  velocity 

-  Tangential  velocity 

-  Dimensionless  tangential  velocity 

-  Axial  velocity 

-  Mean  axial  velocity 

-  Dimensionless  axial  velocity 

*  Transformation  variable 
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X' 


-  First  derivative  of  X  with  respect  to  Z 


Z 

a 

r 

A R 

AX 

F 

v 

v” 

P 

c 

*F 

r*. 

♦ 

Q 


-  Axial  coordinate 

-  Dimensionless  axial  coordinate 

-  Electrostatic  charge  parameter 

-  Circulation 

-  Pennitivity  of  free  space 

-  Radial  change 

-  Axial  change 

-  Viscosity  of  material  constituting  fluid  phase 

-  Kinematic  viscosity  of  fluid 

•  Kinematic  viscosity  of  fluid  phase  material 

-  Fluid  density 

-  Vorticity 

-  Stream  function 

•  Dimensionless  circulation,  vorticity,  and  stream  functions,  respectively 

-  Electrostatic  potential 

-  Dimensionless  electrostatic  potential 

-  Constant  angular  velocity 
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No  of  No  of 

Copies  Organization  Copies  Organization 


2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  HQDA  (SARD-TR) 

WASH  DC  20310-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 


1  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 

Redstone  Arsenal,  AL  35898-5010 

1  Commander 

US  Army  Tank- Automotive  Command 
ATTN:  AMSTA-TSL  (Technical  Library) 
Warren,  Ml  48397-5000 

1  Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSR 

White  Sands  Missile  Range,  NM  88002-5502 


1  Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
Adelphi,  MD  20783-1145 


(ciam.  onty)  i  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 
Fort  Banning,  GA  31905-5660 


2  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-IMI-I 

Picatinny  Arsenal,  NJ  07806-5000 


(Unci***,  only)  i  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 


2  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-TDC 

Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Benet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299-5000 

1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 


1  Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLODL 

Eglin  AFB,  FL  32542-5000 

Aberdeen  Proving  Ground 

2  Dir,  USAMSAA 

ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 
1  Cdr,  USATECOM 

ATTN:  AMSTE-TD 

3  Cdr,  CRDEC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSI 

1  Dir,  VLAMO 

ATTN:  AMSLC-VL-D 
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Intentionally  LEFr  blank. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/uiswers  to  die  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-TR-3207 _ Date  of  Report  February  1991 _ 

2.  Date  Report  Received _ _ 

3.  Does  dds  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ : _ 


- r  .  .  -  -  —  —  -  ■  -  ~~  .  ----- 

'l  ■  1 

5. .  Has  tne  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  Changed  to  improve  future  reports?  (Indicate 
changes td organization,  technical  content,  format,  etc.)  _ _ 


■■■  ■*.  «? - 

A? 

Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


Gty,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 

OLD  Organization 

ADDRESS  _ _ 

Address 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE-— 


Department  of  the  Army 

Director 

US.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21C'<;  -5066 

OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001,  APO,  MD 


POSTAGE  WLL  BE  PAID  8V  ADDRESSEE 


Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


FOLD  HERE- — 


